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ABSTRACT 

Aims. In this paper we present the first on-sky results with the fibered aperture masking instrument FIRST. Its principle relies on the 
combination of spatial filtering and aperture masking using single-mode fibers, a novel technique that is aimed at high dynamic range 
imaging with high angular resolution. 

Methods. The prototype has been tested with the Shane 3-m telescope at Lick Observatory. The entrance pupil is divided into sub- 
pupils feeding single-mode fibers. The flux injection into the fibers is optimized by a segmented mirror. The beams are spectrally 
dispersed and recombined in a non-redundant exit configuration in order to retrieve all contrasts and phases independently. 
Results. The instrument works at visible wavelengths between 600 nm and 760 nm and currently uses nine of the 30 43 cm sub- 
apertures constituting the full pupil. First fringes were obtained on Vega and Deneb. Stable closure phases were measured with 
standard deviations on the order of 1 degree. Closure phase precision can be further improved by addressing some of the remaining 
sources of systematic errors. While the number of fibers used in the experiment was too small to reliably estimate visibility amplitudes, 
we have measured closure amplitudes with a precision of 10 % in the best case. 

Conclusions. These first promising results obtained under real observing conditions validate the concept of the fibered aperture 
masking instrument and open the way for a new type of ground-based instrument working in the visible. The next steps of the 
development will be to improve the stability and the sensitivity of the instrument in order to achieve more accurate closure phase and 
visibility measurements, and to increase the number of sub-pupils to reach full pupil coverage. 

Key words. Instrumentation: high angular resolution - Techniques: interferometric - Planetary systems - Stars: individual: Vega 
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f^ 1. Introduction limit by actively correcting the phas e fluctuations using a de- 
CN formable mirror (iRousset et al.lll990h . Current AO systems are 
^~~^ , The search for faint objects such as exoplanets plays an essential still limited by speckle noise within the central arcsecond 
^ ■ role in the development of high angular resolution techniques. (lGuvonll2005h . and are difficult to implement at optical wave- 
Indeed, the detection of companions orbiting around a parent star lengths. Additional data reduction is necessary to reach high dy- 
requires high resolution on one hand, because of the small sepa- namic ranges at small separations. For instance, the angular dif- 
^ • ration between them, and high dynamic range on the other hand, ferential imaging technique has been app Ued to detect planets 
because of the substantial flux ratio. More precisely, the contrast around HR 8799 dMarois et al.ll2008ll2010l) . and saturated image 
between the companion and the star can easil y reach values from subtracti on has been used for t he detection of a planet around /? 
10'^- 10"'' (for hot Jupiters in the visible, see lDemorv & SeageJ Pictoris dLagrange et al.ll2010l) . Both achieve 4 cr planet-to-star 
1201 lb down to 10"^ or less depending on the wavelength range flux ratios of 10"^ to 5 x 10"^ (in Ks and L' bands) at separations 
and the size and distance of the planet. This kind of target is par- of a few A/D (i.e. a few tenths of an arcsecond). 
ticularly challenging for ground-based telescopes that are sub- In order to reach high dynamic ranges, the AO technique has 
ject to atmospheric effects, the main cause of image degradation, also been coupled with coronagraphic instruments that aim to 
In ideal conditions, without atmospheric turbulence, the angular block the light coming from the star with a mask. Future in- 
resolution of a diffraction-limited telescope is A/D, with A the struments in development should provide contrast ratios on the 
working wavelength and D the pupil diameter of the telescope, order of 10"^, down to 10"^ for the brightest tar gets, at sub- 
Inpresenceof atmospheric turbulence, however, the angular res- arcsecond separation, e.g. SPHERE at the VLT (Beuzit et dj 
olution is ruled by the Fried parameter ro, the equivalent diame- 2010; Mesa et al. 201 1) GPI at Gemini (Macintosh et al. 200^) , 
ter of a telescope working in the diffraction-limited regime. and HiCIAO at Subaru dHodapp et al.ll2008l : ISuzuki et al.ll201( 
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Adaptive Optics (AO) systems can restore the diffraction Impressive results have also been obtained by ISerabvn et al.l 
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(1201 Ol) with one 1.5 m sub-aperture of the Palomar Hale tele- 
scope using a vortex coronograph: they detected the HR 8799d 
planet at a separation of 2/(/D in K, band and announced a 4 cr 
detection limit of a few 10"'' at AjD, very close to the photon 
noise limit. However, this method is difficult to scale to large 
telescopes, where the ratio Djro is much larger 

Post processing of very short exposure images also enables 
the restoration of the diffl^action limit. Speckle interferometry 
(lLabevrielll970h is well suited to detecting binaries but provides 
rest ricted dynarnic ran ge at short separations. Aperture mask- 
ing dHaniff et al.ll987h is another solution exploiting short expo- 
sure images obtained through a mask with holes (sub-apertures) 
placed in the pupil plane in order to remove the noise due to 
atmospheric perturbations. This technique has established itself 
as the standard for high resolu tion imaging with co ntrast ratios 
on the order of 10"^ at AjD dHinklevet alJl201lh . The draw- 
back of this method is that only a small fraction of the tele- 
scope pupil can be used, making photon noise an important fac- 
tor in dynamic range limitation. This can only be overcome for 
bright sources where dynamic range is primarily limited by cal- 
ibration accuracy. Once calibration issues are solved, higher dy- 
namic ranges can be contemplated by increasing the number 
of collected photons. An alternative has been tested at Keck 
Observatory b y dividing the pupil into sets of non-redundant 
configurations (iMonnier et al.l2009h . Differential tilts of the seg- 
mented pupil allow the images of the sets to be focused at dif- 
ferent locations in the image plane. This enables the use of the 
whole pupil along with the aperture masking technique. 

Among the techniques developed to mitig ate atmospheric 
phase fluctuations, the closure phase technique (lJennisonlll958h 
has to be mentioned and highlighted. It consists of combining 
the phases of three different baselines that form a closed triangle, 
and is, by definition, unaffected by phase shifts induced by atmo- 
spheric turbulence. The closure ph ase is a widely use d observ- 
able quantity in apertur e masking (Lacou r et al.ll201 11) and long 
baseline interferometry (IZhao et al. 201 1J. Th is technique seems 
even more promising since Martinachel (|20 1 Ol) has demonstrated 
that in the case of a high Strehl ratio, closure phase-like quan- 
tities called Ker-phases can be measured from direct images. In 
this work, they report on the detection of a known companion 
from Hubble Space Telescope data with a 5 x lO"-' contrast ratio 
at AjD and derive a detection limit of 2 x 10"~ at 0.5 AjD with a 
99 % confidence l evel. 

In this context. [Perrin et all (|2006|) proposed the concept of a 
fibered aperture masking instrument called FIRST for Fibered 
Imager for Single Telescope to achieve high dynamic ranges 
down to a fraction of the diffraction limit. This imager com- 
bines the techniques of aperture masking and spatial filtering us- 
ing single-mode fibers. The fibers only transmit a coherent area 
of the incident wavefront, resulting in a plane wavefront at the 
fiber end. The spatial phase fluctuations are thus converted into 
intensity variations and only residual differential piston terms 
between the sub-pupils remain. Moreover, optical fibers fulfill 
another function: the entrance pupil of the telescope is divided 
into sub-pupils and rearranged into a non-redundant configura- 
tion via the fibers. This non-redundancy is widely used in aper- 
ture masking because all fringe contrasts and phases can be re- 
trieved independently without the attenuation observed when 
complex vectors with random phases are added, as in the case 
of a full pupil. In principle, the combination of single-mode 
fibers and non-redundancy allows the perfect calibration of the 
point spread function over the reconstructed field of view (to 
within photon and detector noise in an exposure time smaller 
than the coherence time of turbulence) and the elimination of 



spec kle noise to reach d ynamic ranges of 10"^ at a fraction of 
A/D dLacour et al.ll2007h . The FIRST concept is therefore com- 
plementary to AO, which allows even higher dynamic ranges but 
at longer distances from the central object. 

A prototype of this instrument wa s developed and tested at 
Paris Observatory (iKotani et al.ll2009l) . Laboratory experiments 
have shown that the image of a simulated binary object could be 
retrieved from the closure phases and visibilities measurements. 
For the sake of simplicity, this prototype uses 9 out of 36 sub- 
pupils, or one quarter of the total pupil area, and works in the 
visible (600nm-800nm). Working in the visible provides good 
resolution (50 mas at 700 nm with a 3-m telescope at the diffi-ac- 
tion limit) and, although the infrared wavelengths are more fa- 
vorable to the detection of low mass companions, could eventu- 
ally offer complementary information by covering a larger spec- 
tral range. 

In this paper, we present the first on-sky results obtained in 
July 2010 with the same prototype set up on the Shane 3-m tele- 
scope at Lick Observatory. In the next section the instrument and 
the key technologies are described. The measurement procedure 
during the observing nights is detailed in Sect. 3 and the data 
reduction method is explained in Sect. 4. Section 5 presents the 
first closure phase and closure amplitude measurements obtained 
on Vega and Deneb. These results and future improvements are 
discussed in Sect. 6. Conclusions are drawn in Sect. 7. 



2. Description of the instrument 

2. 1 . Telescope interface 

The instrument was mounted at the Cassegrain focus of the 
Shane 3-m telescope at Lick Observatory. A special interface 
was constructed to mount the 1.3 m x 0.6 m optical bench, 
consisting of five 50 cm long posts screwed on the bench of 
the AO system and supporting the whole instrument. A beam 
splitter cube, positioned in the light path of the AO system just 
before the wavefront sensor, sends 50 % of the flux from the 
telescope to the FIRST instrument, passing through a 10 cm 
diameter hole in the bench. 



2.2. Optical set-up 

Figure[T]is a schematic representation of the optical set-up. The 
practical implementation is based on a few key components that 
come from advanced technologies: 

- the Segmented Mirror consists of 37 700 jum-diameter 
(vertex-to-vertex) hexagonal segments that can be independently 
positioned in piston/tip/tilt (each segment is controlled by three 
actuators) in order to adjust the phase and steer beams into 
individual fibers. This mirror was manufactured by Iris AO 
(Helmbrecht et al. 2006, 201 1). 

- the Fiber Bundle consists of 36 single-mode Polarization 
Maintaining (PM) fibers (Nufeni PM630-HP, 570 nm cut-off 
wavelength) aligned on a 250yum-pitch hexagonal grid for sam- 
pling of the pupil, without a central fiber because of the tele- 
scope central obscuration. The first prototype uses 9 out of the 
36 available fibers. This bundle was manufactured by Fiberguide 
Industries. 

- the Microlens Arrays consist of 250 jum-diameter fused 
silica microlenses aligned on the same hexagonal grid to fo- 
cus the light into the fibers. They were manufactured by SUSS 
MicroOptics. 

- the v-groove consists of a silicon chip where the PM single- 
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Fig. 1. Schematic view of the optical set-up. F=telescope fo- 
cal point. DP=Dove prism. Ml=pupil mirror. M2=field mir- 
ror PBS=polarizing beam splitter /l/4=quarter wave plate. 
SM=segmented mirror //L=microlens array. FB=fiber bundle. 
VG=v-groove. CYLl, CYL2=cylindrical lenses. P=dispersive 
prism. Inset: side view of the recombination optics. Only three 
fibers are drawn. 
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Fig. 2. a) Entrance pupil with the nine sub-aperture configura- 
tion and the non redundant configuration for recombination. The 
dashed gray circle indicates the central obstruction of the Shane 
telescope. Positions of the fibers in line are (from bottom to top): 
1, 2, 6, 13, 26, 28, 36, 42 and 45. b) Frequency coverage corre- 
sponding to the chosen entrance pupil configuration. 



mode fibers are precisely positioned along a line with a 250 jum- 
pitch, enabling the arrangement of the fibers according to the 
desired non-redundant exit pupil configuration. It was manufac- 
tured by OZ Optics. 

The telescope pupil is divided into sub-pupils by the mi- 
crolens array. Their numerical aperture and hence focal length 
has been chosen so that the injection efficiency into the fibers 
is optimized. The beams corresponding to each sub-aperture are 
focused on the cores of the single-mode fibers gathered in the 
fiber bundle. The segmented mirror is placed in a pupil plane 
and each segment is aligned in order to maximize the injection 
rate into each corresponding fiber An afocal combination of two 



lenses reduces the beam diameters to match the pitch of the fiber 
bundle (700 yum to 250jum). 

The 9-fiber entrance configuration can be seen on Fig. |2] The 
size of the sub-pupils have been chosen in order to maximize 
the sensitivity (which increases with the size of sub-pupils and 
hence decreases with their number if the total pupil size is fixed) 
and optimize the coupling efficiency of the fibers. According to 
ILacou r (2007), a diameter of ~ 3ro is a good compromise. With 
ro on the order of 15 cm at 700 nm under good seeing conditions, 
this leads to sub-pupils of about 45 cm. 

The segmented mirror is used at normal incidence to mini- 
mize flux losses when the segments are tilted. This is made pos- 
sible by a beam splitter which is a polarizing cube used together 
with a quarter wave plate to send half of the flux (one polariza- 
tion state) towards the fibers. The other half is used to control 
the position of the object in the field of view (a CCD camera 
is placed in an image plane). In any case, the two polarization 
states can not be kept because of the birefringence of the fibers. 
Indeed, the polarization states will be phase shifted at the PM 
fiber end, inducing contrast losses. To prevent this, a single po- 
larization state (corresponding to the slow axis of the PM fibers) 
is selected by the beam splitter cube. With this set-up, 100% 
of the light of one polarization is therefore used to produce the 
fringes. 

The fiber outputs are then rearranged according to a one- 
dimensional, non-redundant configuration that can be seen in 
Fig.[2]and is the most compact 9-fiber linear configuration. Each 
pair of sub-pupils forms a unique baseline. A second microlens 
array collimates the 9 output beams. The direction perpendicular 
to the fiber line is dedicated to spectral dispersion. An anamor- 
phic system consisting of the afocal combination of two cylin- 
drical lenses enlarges the beam in the horizontal direction by a 
factor of 20, enabling higher spectral resolution. The beams are 
then spectrally dispersed in the orthogonal direction thanks to an 
equilateral BK7 -prism such that the spectral resolution is around 
150. The final focusing lens forms the image on the EMCCD 
(Electron Multiplying Charge Coupled Device) camera (Luca-S 
EMCCD camera from Andor Technology). The detector con- 
sists of 496x658 lOyum pixels and the quantum efficiency of the 
camera is 50 % at 600 nm, 35 % at 700 nm and 20 % at 800 nm. 
When using the EM mode, the maximum gain is x 200. 

The fibers are therefore the core of the instrument, and also 
the most critical part. Indeed, the total fiber length of each sub- 
aperture must be equal in order to minimize the optical path 
diff'erences (OPD) and thus maximize the fringe contrasts. The 
coherence length gives the maximum OPD where fringes can 
still be observed, and is IQOjjm in air and 67 //m in glass for a 
spectral resolution of 150 at /I = 700 nm. Moreover, equal fiber 
lengths are necessary to minimize chromatic dispersion due to 
the frequenc y-dependent refractive inde x of the fiber and waveg- 
uide eff'ects (JDver & Christensenlll997l) . 

Note that the Dove prism placed in the coUimated beam is 
not absolutely necessary. It provides a way to cover a larger area 
of the (u,v) plane with 9 fibers. Indeed, when the prism rotates 
about the optical axis, the beam rotates at twice the rate of the 
prism. The (u,v) plane coverage is therefore extended by rotating 
the initial frequency positions. However, this component will not 
be needed when the telescope pupil is fully sampled. 

3. Principle of measurements 

The critical parts of the alignment are the injection into the 
single-mode fibers and re-collimation at the v-groove fiber ends. 
The core diameter of the fibers is 4 yum, while the focal length of 



E. Huby et al.: FIRST, a fibered aperture masking instrument 



the microlens is about 1 mm. Stable alignments and a systematic 
procedure during night observations were thus needed. 



3.1. Injection optimization procedure 

The segmented mirror is the key component for fiber coupling 
optimization. Software has been developed to control each seg- 
ment. A raster scan has been implemented to automatically find 
the best orientation of each segment by measuring the flux trans- 
mitted by each fiber onto the EMCCD camera as a function of 
segment tip-tilt values. This alignment is aimed at correcting the 
static aberrations introduced by the optical train, while the IR- 
optimized AO system actively corrects the turbulent fluctuations 
of the wavefront. At the scale of each sub-pupil, the correction 
amounts to removing part of tip-tilt while higher order aber- 
rations are left unchanged. In addition, the average differential 
phase from one pupil to another remains negligible as the AO 
system reconstructs a continuous wavefront, so the AO system 
is equivalent to a multi-pupil tip-tilt and differential piston cor- 
rector In addition, a global image alignment was performed with 
a motorized steering mirror (Ml on Fig. [1). Pointing errors are 
monitored on the CCD placed in an image plane. The alignment 
is adjusted by maintaining the target at a fixed position in the 
field of view. 

Details of the instrument throughput are presented in Table[T] 
The total efficiency is deduced from data taken on Vega by com- 
paring the number of photons expected from the stellar spec- 
trum in the 600- 760 nm range to the number of counts mea- 
sured on the detector in 200 ms. Telescope, AO transmission and 
FIRST optics transmission have been estimated with the theo- 
retical transmission and reflection coefficients. An upper value 
of the fiber injection efficiency has been obtained by assuming 
that injection losses are du e to residual tip-tilt, mode mismatch 
(a 22% loss according to IShaklan & Roddieij[T988l) . and high 
order phase perturbations. 

Tip-tilt data were regularly measured by the AO wavefront 
sensor during the observations to evaluate the typical perfor- 
mance of the tip-tilt correction. The wavefront sensor consists of 
sub-pupils of the same diameter as the FIRST sub-pupils (that 
is, about 43 cm on the primary mirror). The standard deviation 
of the residual fluctuations is on the order of 125 mas per sub- 
pupil on the sky, or a l.Ojum lateral offset of the Airy pattern 
on the fiber core. The average injection efficiency due to tip-tilt 
fluctuations is of 76 % at 700 nm, computed with the formula: 



;; = exp 



-20V' 



(1) 



with 6 the tip-tilt angle on the microlens, / the microlens focal 
length, and cjq the waist of the fiber mode. This expression can 
be used under the assumption of perfect mode matching between 
the beam and the fiber mode. 

At the time of the observations, the typical ro at 550 nm was 
22 cm, and hence the spatial phase variance over a sub-pupil 
was 1.9rad" at 700 nm. The phase variance (free of tip-tilt) is 
0.25 rad^ which translates to a coherent energy of 78 %. Taking 
into account the mode mismatch and residual tip-tilt, the fiber 
injection efficiency is 46 %. 

The throughput budget yields a theoretical efficiency of 
1.2%, which overestimates the measured efficiency of 0.21 % 
by a factor of ~5.7. This difference may be explained by im- 
perfect alignments and also non-ideal performance of the optics 
(ageing, dust, etc ... ). For instance, assuming that all 46 sur- 



Table 1. Instrument transmission eflSciency. 

Item Transmission ( %) 

Telescope and Adaptive Optics 48 

FIRST optics 22 

Entrance beam splitter 50 

Polarization selection 50 

Injection rate into fibers 46 



Theoretical efficiency 
Measured efficiency 



1.2 
0.2 



faces (that are exposed to air) share the same contribution to the 
total loss and applying an additional efficiency of 97 % per op- 
tical surface makes the theoretical transmission drop to 0.25 %. 
This may explain the poor efficiency of FIRST during the first 
tests reported in this paper This clearly leaves some margin for 
improvement. It is obvious that the current sensitivity limit can 
potentially be improved by reducing the number of optical sur- 
faces thanks to a different set-up and by improving the quality of 
optics. Achieving the theoretical 1 % overall efficiency therefore 
seems a sensible goal and would lead to a sensitivity limit im- 
proved by 1 magnitude in the photon noise-limited regime. This 
is further discussed in Sect. |6] 

Note that the data reported here were obtained with only four 
fibers at sufficiently high signal-to-noise ratios to produce stable 
fringe patterns (thus leading to stable and relatively accurate clo- 
sure phase measurements, see Sect. 4). Two out of the nine fibers 
fell inside the central obstruction of the telescope, as shown in 
Fig |2] (this will be changed in the future with a new fiber con- 
figuration). In addition, one fiber was damaged during transport 
and did not transmit light. Two other fibers transmitted only low 
fluxes or produced low contrast fringes (due for example to a 
high flux ratio between some pairs of sub-apertures, or perhaps 
to unequal fiber lengths causing the mean OPD to be comparable 
to the coherence length). 

3.2. Image acquisition 

After alignment of the nine segments, data cubes of 50 images 
were acquired, with integration times of 200 ms. Although 
optimized for the infrared, the AO system provides sufficient 
wavefront corrections to stabilize the phase of the fringes, thus 
fulfilling the role of a fringe tracker and making longer inte- 
gration times possible (coherence time in the visible is only on 
the order of 10 ms). The choice of integration time is the result 
of a compromise between sensitivity and fringe contrast loss 
that has been roughly evaluated during observations. The study 
of the residual piston between sub-pupils after the correction 
applied by the AO system is not trivial and needs to be detailed 
in a longer study that is beyond the scope of this paper 

The camera was used in the Electron Multiplying (EM) 
mode with the highest gain (x200), making the readout noise 
negligible with respect to the other noise sources (the dark 
current and the photon noise are amplified, while the readout 
noise is not). However, the noise is amplified by a factor of 
F - V2 when using the EMCCD mode. This is equivalent 
to a reduction of the si gnal-to-noise ratio by a factor of yl 
dDenvir & Conrovll2003h . while the average value of the signal 
itself is not affected. Since the dark current is insignificant in the 
visible, the acquisitions are photon noise-limited, apart from the 
multiplicative noise factor 

A full acquisition was completed in six steps: four sets of 
images on the object were taken for the different positions of 
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Table 2. Observation log. 



Target 


Start time (UT) 


Acquisition 


Dove position (°) 


Vega 


06:01:17 


10 datacubes 


45 


Vega 


06:09:16 


10 datacubes 


45 


Vega 


06:25:12 


10 datacubes 


90 


Vega 


06:32:06 


10 datacubes 


90 


Vega 


06:43:20 


10 datacubes 


90 


Vega 


06:51:23 


10 datacubes 


90 


Vega 


07:09:11 


10 datacubes 





Vega 


07:17:00 


10 datacubes 





Vega 


07:33:48 


10 datacubes 


135 


Vega 


07:41:31 


10 datacubes 


135 


Vega 


07:51:21 


calibration files 




Vega 


07:58:21 


sky background 




Deneb 


08:22:20 


10 datacubes 


135 


Deneb 


08:31:20 


10 datacubes 


135 


Deneb 


08:51:16 


10 datacubes 


90 


Deneb 


08:58:35 


10 datacubes 


90 


Deneb 


09:06:51 


10 datacubes 


90 


Deneb 


09:28:05 


10 datacubes 


45 


Deneb 


09:37:44 


10 datacubes 





Deneb 


09:44:58 


sky background 





the Dove prism (0°, 45°, 90° and 135°), followed by one set 
of sky background with telescope offset 30arcsecond and a 
calibration sequence. During this sequence, the flux of each 
sub-pupil was independently recorded, with all segments but one 
being deliberately misaligned. This step is aimed at measuring 
individual fiber images, which are necessary for modeling the 
interferogram (see Sect. 4). 

For the first on-sky tests, the target list was limited to bright 
objects with Rmag ^ 3 near zenith as the instrument suffered 
from mechanical flexure. In this work, we present the results 
of observations of Vega (AOV star, Vmag^O.OB, Rmag=0.1) and 
Deneb (A2I star, Vmag=l-25, Rmag=l-14) as a calibrator. The 
observations were conducted on July 30 2010 from 06:01 UT 
to 09:45 UT under good seeing conditions, with a visible ro 
estimated to be 22 cm at 550 nm on average by the Lick-AO 
system. The observation log is reported in Table |2] indicating 
the number of data cubes (each containing 50 frames) and the 
position of the Dove prism. 




100 100 

Spectral Channel [pixel] 



-100 100 

Group Delay [fim] 



Fig. 3. The three images to the left are samples of the acquired 
data on Vega, where fringes are particularly visible. The right 
panel is the 2D Spectral Power Density averaged over 40 data 
cubes of 50 images each, in a logarithmic scale. One can distin- 
guish 13 peaks, meaning that at least 6 fibers contributed to the 
interferograms in this particular set of data. 



iLacour et al.l (l2007h . the observables are the complex contrasts 
ju,y for each pair of sub-pupils / and j: 



idj-^i+ifij) 



(2) 



//,V = VijGiG* = |y,7l|G;||G;|e' 

with Vij the complex visibility (modulus \Vij\ and phase ifij ) 
and Gi the complex transmission in the pupil ; (modulus |G,|, the 
transmitted flux, and phase O,). The ///y are basically the mea- 
sured Fourier components, while the Vij are the intrinsic Fourier 
components of the object. The phases of the complex transmis- 
sions Gi account for residual differential piston (including atmo- 
spheric and instrumental OPD). 



4. 1 . Model for the interferograms 

The acquired data are spectrally dispersed. Each column of the 
image corresponds to one spectral channel of width 6A. The in- 
tensity recorded for the channel at wavenumber cr is: 



1=1 



4. Data reduction ^r 

Typical images acquired with FIRST are presented in Fig.|3] The ^ W - 2_i '^'' '^i^-"-' 
wavelengths are in abscissa and the OPD is in ordinate (pre- 
sented as the position x on the detector). Each image consists of 
the 133 spectral channels used for the data reduction. They cor- 
respond to the 399 pixels with the highest signal-to-noise ratios 
in the raw images, binned in sets of 3 to increase the signal-to- 
noise ratio. The 2D spectral power density shows that the images 
are the superposition of at least 13 fringe patterns. However, this 
is not part of the data reduction that is presented in this work, 
since the data are processed directly in the image plane rather 
than in the Fourier plan e, as for the AMBE R long-baseline in- 
terferometry instrument (iMillour et al. 1120041) . 

Indeed, because of the multi-axial beam combination scheme 
with spatial coding and spectral dispersion, the final interfer- 
ogram is a linear combination of fixed patterns with parame- 
ters linked to injected fluxes, fringe contrasts and phases. If the 
carrying waves are well determined, the fringe parameters can 
then be measured through model fitting. Using the formalism of 



2 J ^ai(x)aj(x)%{fiije^""'''"i), 

></<A'p„p 



(3) 



where x is the position on the detector (x takes discrete values 
xii, multiples of the pixel size of lOjum), A^pup is the total number 
of sub-pupils, «B the number of baselines, and D the distance to 
the image plane. ai(x) is the normalized envelope function for 

np 

the i pupil such that 2 ^iixk) - |G,|, «p being the number of 

k=l 

pixels, and B,y is the baseline formed by the pupils / and j in the 
exit pupil. 

This expression can also be written: 

I{x) = Yj\Gi\Wx) 
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I 

2 2_] ^ (/"ij) ^ai{x)aj{x) cos [infijX^ 



j<i<N„ 



I 

-2^3 [/Jij) ^ai(x)aj(x) sin [ijrfijx) , (4) 

j<i<N„ 

with fij = cr-^, the spatial frequency corresponding to the base- 
hne ij. Assuming the photometric fluxes |G,| and the enve- 
lope functions a,(x) are known, the interferogram can be cor- 
rected from the fixed component and becomes: 



Rij and 7,^ are the best fit parameter sets in the least-squares 
sense. 

This simplified matrix inversion can be written under the 
assumptions of uncorrelated pixels and constant variance over 
all pixels. In practice, the matrix M is rectangular and the 
generalized inverse matrix is determined via singular value 
decomposition . Acc ording to the formalism introduced by 
Millour et al.i (l2004h . the matrix M containing the carrying 
waves is called the visibility to pixel matrix (V2PM), while 
the inverse of this matrix is called the pixel to visibility matrix 
(P2VM). 



-■pup 

J'(x) = I{x) - J] \Gi\'a'(x) 
1=1 

«B «B 

= 2 ^{^'u)ciM)- 2 3(juo>,7(x), (5) 



with 



Cij{x) — 2 J ai(x)aj(x) cos yZTrfijx] , 
Sijix) = 2 ^ai(x)aj{x) sin (in fijxY 



(6) 
(7) 



Cij and Sij are defined to be the carrying waves that constitute the 
model of the fringe pattern. They have to be precisely calibrated 
in order to retrieve the real and imaginary parts of the /i,j by 
fitting the fringes with this model. 

4.2. Calibration 

The envelope functions of each fiber, fl,(jc), are defined from the 
calibration data described in the previous section. The spatial 
frequencies, that depend on the wavenumber o- and the base- 
line Bij, also have to be determined. Assuming that the base- 
lines are known (the v-groove pitch is given with a 0.5 /im pre- 
cision), a first approximation of the central wavelength of each 
channel is made using one spectral absorption line observed in 
the stellar spectrum (//„ at 656 nm) and two lines due to the at- 
mosphere (O2 at 687 nm and 760 nm). More accurate values of 
the wavenumbers and baselines are obtained by fitting interfer- 
ograms in each spectral channel. The resulting wavelengths are 
finally fitted with a 5* order polynomial function over the spec- 
tral channels. 



4.3. Fringe fitting 

Once this calibration is done, the problem is linear regarding the 
frequency components and can be expressed in terms of matri- 
ces. For one spectral channel, the interferogram is sampled by 
the number of pixels rip: 



( I{xi) 



•. liXtif) J 
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fRij^ 
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(8) 



with M a «p X 2 «b matrix defined by the carrying waves: 

n 



M = 



'1 



,1 



"B 1 

'-np ■'«p 



'1 



^«o / 



(9) 



4.4. Measured quantities 

Fringe contrasts yUy are measured with: 



(10) 



According to Eq. |2] the individual fluxes |G,| are needed to re- 
trieve the object visibility, often measured through the squared 
visibility: 



Wij 



lA'oi- 



(11) 



|G,'P|G/- 
Closure phase is by definition the phase of the bispectrum 



f^ijk = (MijMjkMik*)' 



(12) 



with ijk referring to the pupils forming the considered baseline 
triangle and () indicating the average over all measurements. The 
most interesting property of this quantity is that it is indepen- 
dent of differential pistons, which naturally cancel each other out 
when summing the phases of three baselines that form a closed 
triangle: 



CPijk = arginiji) = ipij + (fjk - ipik, 



(13) 



with (fij the phase of the complex visibility V,y. 

Since the squared visibility and closure phase are expressed 
as a power law of a noisy quantity (additive photon noise), unde- 
sired photon biases affect these estimates. High pre cision mea- 
surem ents require bias-free estima tors (see Goodman & Belsherl 
(1 19761) for squared visibility and IWirnitzen (Il985l ) for bispec- 
trum.) Such estimators have not been used for the reduction of 
the first data. Indeed, only closure phases of unresolved targets 
have been retrieved for now and in this case photon biases cannot 
be the cause of non-zero closure phases (the bias is necessarily a 
positive real number). Concerning the visibilities, Eq. [TT] states 
that the individual fluxes of each sub-pupil are needed to esti- 
mate the visibility terms. However, the FIRST instrument does 
not include any photom etric channel to measure these quanti- 
ties. iLacour et al.i (l2007h developed an algorithm to retrieve ob- 
ject visibilities directly from the fijj measurements. This is pos- 
sible if the rank of the matrix linking the observables (//,j) and 
the unknowns (G, and Vtj) is at least equal to the number of un- 
knowns, that is, if the entrance pupil configuration shows a suf- 
ficient number of redundancies (which is the case for the chosen 
pupil configuration illustrated in Fig.|2). 

For the first-light observations, only four fibers could be 
used to measure closure phases with reasonable precisions (see 
Sect. 3.1), reducing the pupil to a non-redundant configuration 
and making the problem ill-posed (12 complex equations for 
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Table 3. Mean closure phase results and statistical errors. 
Baseline Triangle CP Vega (°) CP Deneb (°) Calibrated CP (°) 



1: 1-2-6 


-10.4 + 0.4 


-9.4 ±0.4 


-1.0±0.6 


2: 1-2-13 


3.5 ±0.3 


1.5 ±0.4 


2.0 ± 0.5 


3: 1-6-13 


25.5 ± 0.4 


23.6 ±1.1 


1.9 ±1.2 


4: 2-6-13 


12.2 ±0.4 


10.9 ±0.6 


1.3 ±0.7 



Table 4. Comparison of statistical errors. The mean error bars 
correspond to the average of all individual error bars (shown for 
each measurement point in Fig.|4]i. The mean dispersions are the 
standard deviations of all measurement points. 



Vega 



Deneb 







Mean 




Mean 


Baseline 


Mean 


dispersion 


Mean 


dispersion 


triangle 


error bar (°) 


of data (°) 


error bar (°) 


of data (°) 


1: 1-2-6 


3.5 


3.5 


3.0 


3.0 


2: 1-2-13 


3.3 


3.9 


5.4 


3.5 


3: 1-6-13 


4.0 


3.9 


7.3 


5.0 


4: 2-6-13 


4.0 


3.8 


4.8 


3.9 



6 complex visibility terms and 4 complex transmission terms). 
Calibrated fringe contrasts could therefore not be measured from 
the first run data. 

Another quantity can be built from the uncalibrated contrasts 
in the same spirit as closure phase: 



CA 



ijki 



lj"i./lly"»l 
It^ikWl^jll 



VijVkl 



V„Vj, 



(14) 



This quantity is called closure a mpUtude (CA) and onl y depends 
on the moduli of the visibilities (iReadhead et alJfl980h . 



5. First results 

5.1. Closure Phase results 

The first closure phases are calculated on Vega using Deneb as a 
calibrator (Table |3). Both are supposed to be unresolved targets 
with achromatic characteristics. The closure phases are defined 
by the fibers that form the corresponding triangle. As mentioned 
above, only four fibers have been used to compute these four 
closure phases: 1, 2, 6 and 13 are their positions in the v-groove 
(see Fig.|2|i. The longest baseline on the primary mirror is 2.3 m 
and corresponds to the pair 6-13. 

The raw closure phases are presented in Fig. 2] Each value 
is computed from the mean of 50 images (10 s of observation, 
the exposure time for each image being 200 ms). Bispectra are 
estimated for each of the 133 spectral channels and averaged 
to yield wideband bispectrum values {AA =160 nm centered on 
Aq - 680 nm). Raw statistical error bars for each set of 50 images 
are first derived from the standard deviation over the 50 meas ure- 
ments weighted by the modulus of the bispectrum. Second, a 
constant value (the average) is fitted to the sequence of raw mea- 
surements and a reduced x~ is calculated. If statistical error bars 
underestimate the fluctuations of the data then the x~ is larger 
than 1 . In this case, individual variances are rescaled through a 
multiplication by the x^ value (to get;^'^=l) to derive a more 
realistic error bar for the closure phase estimates of Table [3] 

The modulus of the bispectrum, |yU,yt|, is an effective indicator 
for the detection of the fringes. A signal-to-dark ratio (SDR) has 
been defined by: 



SDRijk = 



iMijkl .iq3.^ 



If^ijkdarkl 



(15) 



with luijkdarkl the mean modulus of the bispectrum obtained by 
applying the pixel to visibility matrix to sky background images, 
and Rmug the magnitude of the target. The multiplying factor in 
the second part of the equation allows for comparable SDR lev- 
els for targets with diff'erent magnitudes. This quantity does not 
exactly define a signal-to-noise ratio but constitutes a way to dis- 
criminate the data. The data with SDR < 210 have been rejected; 
this threshold has been chosen in order to keep around 30% of 
the data in the worst case (i.e. for triangle 1-6-13 on Deneb). 
The holes in the sequence of data are due to this selection, espe- 
cially for Deneb. Moreover, there are gaps with no measurement 
points that are represented by gray rectangles: they correspond 
to the optimization of the flux injection into the fibers, requiring 
about 30 seconds per segment of the mirror (see Sect. 13.11 . 

Datacubes were taken for different positions of the Dove 
prism (see Table|2]i. Rescaled error bars in Fig.2]are comparable 
to the standard deviation of all measurement points, as shown 
in Table |4] This means that there is no significant drift or sys- 
tematic error on timescales of 1-2 hours. In other words, the 
closure phase measurement is unaffected by mechanical flexures 
or Dove prism rotation. 

We have decided to average results for different positions 
to improve the SDR. The SDR for individual positions was not 
high enough to detect noticeable effects and we have assumed 
that 1) instrumental biases would be identical and 2) that any as- 
trophysical signal would be small. The raw closure phases reach 
high values (up to 26°) and are of instrumental origin. The bis- 
pectrum biases due to photon and detector noise could not be 
properly subtracted because of the lack of calibration data se- 
cured during this first run. This is why we have selected only 
high SDR data to reduce the impact of noise biases. In any case, 
such biases cannot be the cause of non-zero closure phases. 

Final calibrated closure phases are obtained by subtracting 
the Deneb from the Vega closure phases. The results are pre- 
sented in Fig. |5] and listed in Table [3] The statistical error is 
on the order of 1° and at best 0.5°. However, non-zero closure 
phases are measured at the 2 to 3-cr level. It is hard to say at 
this stage if these are uncalibrated systematic errors or of some 
astrophysical origin. Vega is known to be surrounded by a de- 
bris disk, which may contribute to visibilities in scattered light. 
Systematic er rors on closure phas es are also encountered in aper- 
ture masking (iLacour et al.ll201 u) at the 0.3° level and still need 
to be understood. In our case, they could be due to polarization 
errors (misalignment of fiber neutral axes, for example) and need 
to be investigated. 



5.2. Closure Amplitude results 

As explained in Sect. 4.4, visibility amplitudes could not be ac- 
curately measured since beam photometries cannot be retrieved 
without redundancy in the entrance pupil. An order of magni- 
tude can however be estimated assuming typical flux values for 
the individual beams. The mean total flux per spectral channel 
is approximately 2x10^ ADU, leading to a typical flux value of 
5x10"* ADU per fiber With a mean coherent flux ranging be- 
tween 1x10^ and 4x10^* ADU, we estimate the raw contrasts for 
the 6 baselines to roughly range between 10 % and 40 %. 

Closure amplitudes (CA) have been measured and are pre- 
sented in Fig.|6] These estimates have been obtained by follow- 
ing the same procedure as for the closure phase measurements : 
computing estimates by selecting the data according to a given 
threshold, estimating the mean over all points, and correcting 
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Fig. 4. Estimates of the closure phase of Vega (blue squares) and Deneb (red triangles) with rescaled error bars, corresponding to the 
contributions of four sub-pupils numbered 1, 2, 6 and 13: a): 1-2-6 ; b): 1-2-13 ; c): 1-6-13 ; d): 2-6-13. Each estimate corresponds 
to 10 s of observation (average of 50 images of 200 ms integrations). The error bars represent the standard deviation of the mean 
over the 50 frames (i.e. the frame-to-frame standard deviation weighted by the modulus of the bispectrum divided by VSO). The 
gray rectangles represent the periods when no data were taken. 
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Fig. 5. Top: Closure phase results for Vega (blue squares) and 
Deneb (red triangles), which are both supposed to be unre- 
solved targets. The error bars are comparable to or smaller 
than the symbol size. The exact values are given in Table |3] 
Bottom: Calibrated closure phase of Vega. The error bars are the 
quadratic sum of errors from Vega and Deneb. 



Fig. |6] shows that the closure amplitudes are stable over about 
half an hour. The observation procedure could therefore be op- 
timized by alternating the target and the calibrator observations 
more frequently. 

The calibrated closure amplitudes of Vega are shown in Fig. 
|7] They are close to the expected value of unity for these un- 
resolved targets. CAi and CA2 are measured with precisions of 
about 25 %, which is not sufficient for high dynamic range de- 
tection. According to Eq. [16] these closure amplitudes have 2 
baselines in common: (2-1) and (13-6). The biases observed in 
the resulting closure amplitudes could possibly come from the 
contrast measurements of these bases, but as for closure phases, 
we must further investigate the sources of these biases. 

On the other hand, CA3 is stable: dispersion of the data 
points in Fig. |6]is 0.2 whereas the overall mean error bar is 0.4. 
Moreover, the estimates of the mean are 1 .0 for Vega and 1 . 1 for 
Deneb, with accuracies of 10 %. In other words, the third closure 
amplitude is fully consistent with a point source nature for both 
objects and a proper data analysis. 



5.3. Dynamic range 



the final error by the x^ calculation. The closure amplitudes are 
numbered from 1 to 3: 



CAi = 



/^a-i/^B-e 



^6-1^13-2 



,CA2 = 



H2-mn-6 



^^\i-\^^(,-2 



, CA3 = 



P-6-\P-n-2 



l^l3-\f^6-2 



.(16) 



Given the results presented in the previous Section, only the clo- 
sure phases will be used to evaluate the dynamic range detection 

limits. 

iLacour et al.l ( 1201 lb gives a 1 cr dynamic range estimation for 
sparse aperture masking using 7 sub-pupils, as a function of the 
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Fig. 6. Estimates of the closure amplitudes of Vega (blue squares) and Deneb (red triangles) with rescaled error bars, corresponding 
to the contributions of four sub-pupils numbered 1, 2, 6 and 13: a): CAi ; b): CA2 ; c): CA3. Each estimate corresponds to 10s of 
observation (average of 50 images of 200 ms integrations). The error bars represent the standard deviation of the mean over the 50 
frames (i.e. the frame-to-frame standard deviation divided by VSO). The gray rectangles represent the periods when no data were 
taken. 



error on closure phase measurements based on a Monte Carlo 
simulation: 



Icr dynamic range = 



400 



(r(CP)A 



(17) 



This equation can be extrapolated to estimate the Act dynamic 
range for any number of baselin e s wb ( based on the theoretical 
derivation by [Baldwin & Haniffl (l2002h . who have established 
that the dynamic range increases with the square root of the num- 
ber of baselines): 



4cr dynamic range ~ 22 



V«B 



cr{CP). 



(18) 



degr 



The calibrated closure phase estimations presented in this 
first paper reach a final accuracy of 0.75 ° on average, resulting 
from 17 min total integration time on Vega and 12 min on Deneb. 
This leads to a Act dynamic range of 180 in the reconstructed 
image, if all 9 fibers were equally transmissive. 

5.4. Sensitivity 

During the first run, the target list was limited to bright objects. 
Nevertheless, the sensitivity of the instrument can be estimated 
by extrapolating the results obtained on the faintest target, which 
is Alf And (B9II star, V,nag=2.06, Rniag=2.1), observed under the RmagUm = RmagO + 2.5 log 
same conditions (200 ms integration time). The same procedure 
has been applied to the 4 sets of 10 datacubes acquired on this 
target to obtain closure phase measurements. The statistical error 
bar of one acquisition (resulting from 50 200ms-images) is 6.6 ° 
on average, and 4.4 ° for the triangle with the highest signal-to- 
noise ratio. When combining all data, the final statistical error is 
0.6 ° for 7 min effective integration time. 

Sensitivity limits can be extrapolated from the results ob- 
tained on the different targets and can be derived for a given 




1 2 3 

Closure Amplitude Number 

Fig. 7. Top: Closure amplitude results for Vega (blue squares) 
and Deneb (red triangles). Boffom.Calibrated closure amplitude 
of Vega. The error bars are the quadratic sum of errors from Vega 
and Deneb. 



closure phase accuracy, measured on targets with R,„ago and in- 
tegration time tq: 



To 



(19) 



These limits are derived under the assumption of the photon 
noise limited regime. A signal-to-noise ratio of at least 1 is 
needed per exposure and per spectral channel in order for the 
overall integrated signal-to-noise ratio to be reasonable. In other 
words, at least one photon is needed per fiber per frame per 
coherence time. As FIRST is working behind an AO system, 
the coherence time is theoretically infinite and integration times 
could be set to seconds or more for the faintest targets. 
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Fig. 8. Sensitivity limits as a function of the integration time 
for given closure phase accuracies per acquisition. These curves 
have been extrapolated from the accuracies measured on Vega 
(square at R„,flg=0.1), Deneb (triangle at R„,fl^=1.14) and Alf 
And (cross at R,„ag=2.l). Two curves are drawn per target: the 
black one is extrapolated from the mean accuracy per acquisi- 
tion and the gray one from the final best statistical error The 
accuracies are given for uncalibrated measurements. 



Sensitivity limits are plotted in Fig. [8] as a function of ef- 
fective integration time for one acquisition (SOxtq) and a given 
accuracy. These accuracies are for uncalibrated measurements, 
meaning that the final accuracies should be multiplied by a fac- 
tor y/2 in the best case (reference target of the same magnitude as 
the science target), and effective integration times should be dou- 
bled to take the reference target observations into account. The 
current state of the prototype could measure closure phases of 
targets at Rmag ^ 4 with a 3.7 ° accuracy per acquisition (uncal- 
ibrated) in about 6 minutes, and with a 0.6 ° accuracy in 40min. 
The achievable magnitude is higher for relaxed constraints on 
the accuracy. 



6. Discussion 

The first on-sky results with FIRST demonstrate that it is possi- 
ble to operate this type of instrument with multiple single-mode 
fibers. Several key ingredients have been used: the segmented 
mirror that precisely steers beams into fibers, and the microlens 
arrays, and fiber bundle which compactly focus light into the 
fibers. However, the throughput is low and the accuracy of the 
closure phase measurements are currently limited by system- 
atic errors. Future developments of the FIRST instrument will 
be dedicated to improving these aspects. 

6.1. Sensitivity and robustness: Long-term improvements 

Given the throughput budget presented in Table \T] some im- 
provement can be made by decreasing the number of optical 
surfaces to reach a higher theoretical efficiency (currently a few 
percent). 

This could be achieved in several ways. First of all the optical 
design could be simplified insofar as it could be better integrated 
into the AO system. Indeed the prototype set-up is complex be- 
cause the AO and FIRST are separate systems. For instance, the 



mirrors required to control the pupil and image planes in FIRST 
(Ml and M2) are redundant with their counterparts in the AO 
system. This is the same for the optical train required to extract 
the beam from the AO system (one beam splitter cube and two 
mirrors, which are not drawn in Fig. [1] are currently needed). 
The set-up could also be simplified by placing the segmented 
mirror at quasi normal incidence. This would cause negligible 
flux losses compared to the gain in throughput: the polarizing 
beam splitter cube and the quarter wave plate could be removed 
(removing 10 reflections/transmissions). In this case, the control 
of the image in the field of view would not be done continuously 
but alternately, by putting miiTorMl or M2 on a translation stage 
to let the beam reach the alignment CCD. 

Second, part of the output optics of FIRST could be re- 
placed by a 3D integra ted photonics pupil remapper (as for the 
Dragonfly instrument, i Tuthill et al.l gOly), avoiding the fiber- 
interfacing between the fiber bundle and the v-groove, and pro- 
viding a more robust and thermally stable instrument. Integrated 
optics could also replace some of the recombination optics, such 
as the anamorphic system (see 'Benis tv et al.l ( l2009l) for ID in- 
tegrated optics and Minardi & Pertschi (l2010h for 2D integrated 
optics). The development of these technologies is in progress. 
Currently only prototypes working in the near-infrared wave- 
length range have been tested on-sky, while integrated optics at 
visible wavelengths remain to be developed. 

6.2. Next prototype: FIRST-1 8 

In order to reach better performance and allow for observations 
of astrophysical interest in the shorter term, three axes of im- 
provement must be implemented on the current prototype. A 
major identified limitation to the measurement accuracy is me- 
chanical stability. The prototype was first made to be a labo- 
ratory instrument, yet was mounted directly on the Cassegrain 
focus of the Shane telescope. Mechanical instabilities such as 
flexure have turned out to be a critical issue during the first run. 
The bench will undergo substantial improvements to provide 
better injection into fibers and therefore increase the sensitiv- 
ity. Second, the next version of the prototype, named FIRST- 18, 
will consist of 18 fibers, thus increasing the number of available 
baselines. The pupil mapping will be updated to avoid the central 
obstruction of the telescope. Finally, better protection from dust 
and a better coating of the optics will ensure a higher throughput 
to reach the 1 % theoretical efficiency. 

As stated by Eq. [18] the dynamic range will be increased 
due to the larger number of fibers and improved sensitivity and 
mechanical stability, enabling more accurate closure phase mea- 
surements. The cuiTent accuracy on closure phases may be pri- 
marily limited by systematics which will be reduced by data pro- 
cessing and by better calibration procedures. 

7. Conclusion 

We have obtained the first on-sky results with FIRST, a fibered 
aperture masking instrument that aims at providing a unique 
combination of high contrast and high angular resolution. 
Although the described performances are not very high, better 
precisions and sensitivity are expected in a near future. 

Apart from the systematic errors that may remain, we have 
measured calibrated closure phases on Vega with an accuracy 
of 0.5 ° in the best case by co-adding ~ 100 acquisitions (with 
accuracies around 3-4° per acquisition of 10 s), totaling 17min 
exposure time. This leads to a 4cr dynamic range of around 260 
in a reconstructed image when considering a 9-fiber prototype. 
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This dynamic range w ould increase with the square root of the 
number of acquisitions (iPerrin et al.ll2006l) . if not Umited by sys- 
tematic errors. By extrapolation, the 4cr dynamic range could 
therefore be on the order of 10^ when co-adding 14 IVmin- 
acquisitions, equaling about 4 hours total effective integration 
time for both target and calibrator However if we achieve a clo- 
sure phase accuracy of 0.4 ° in 17 min, the same dynamic range 
could be achieved in 2.6 hours integration time. Moreover, when 
the 30-fiber instrument will be completed, the number of base- 
lines will be 210 (if we assume 2 recombination paths of 15 
fibers), and a dynamic range of 10"* could be reached in less than 
3 hours if a bias-corrected closure phase accuracy of 0. 1 ° is 
achieved in 17 min. 

Closure amplitudes have not been used to evaluate the dy- 
namic range detection limits because of the limited precision on 
these estimates (10 % in the best case). With a fully working in- 
strument, data processing would allow the computation of the 
visibilities, leading to better performance in high dynamic range 
detection. 

The preliminary data presented in this paper, although they 
do not show very high performances, constitute an important 
demonstration of the fibered pupil remapping concept applied to 
a monolithic telescope. FIRST in the visible range is not compet- 
itive resolution-wise with long baseline interferometers. Long- 
baseline instruments provide much higher angular resolution in 
the visible, but with more limited field of view, unless spectral 
resolutions of typically a few thousands are used. FIRST on 3 to 
10-m class telescopes is a good complement, as it gives access 
to a range of large objects and to spatial frequencies otherwise 
unaccessible with either IR-optimized adaptive optics systems 
or long baseline interferometers. Moreover, a fibered aperture 
masking instrument is also complementary in terms of dynamic 
range since it can provide high dynamic range at the diffraction 
limit, while AO systems can reach high dynamic range at a few 
AjD. Possible science applications are stellar surfaces for giants 
and supergiants, circumstellar environments, debris disks, etc ... 
Although FIRST could eventually be used to observe the clos- 
est and brightest planets, a near-infrared version would be better 
suited to exoplanetary system studies. 

Our future efforts in the development of FIRST will be con- 
centrated on the improvement of the sensitivity along with the 
accuracy of the closure phases and visibilities. Dynamic ranges 
up to lO-'-lO^ and a sensitivity limit of R,„„p > 4 seem within 
range in the not too distant future. 
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